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Abstract 

A method to estimate phase velocity and attenuation of acoustic waves in 
the presence of liquid water in a snowpack is presented. The method is based 
on Biot’s theory of wave propagation in porous materials. Empirical relations 
and a priori information is used to characterize snow as a porous material as 
a function of porosity. Plane wave theory and an equivalent pore fluid are 
used to solve Biot’s differential equations and to asses the impact of the air 
and water in the pore space. The liquid water in the pore space of a snow 
pack reduces the velocity of the first compressional wave by roughly 300 m/s 
for every 0.1 increase in liquid water saturation. Also the attenuation of the 
compressional waves is increased with increasing liquid water content. Two 
end member models for compaction are evaluated to asses the importance of an 
independent density measurement for an estimate of liquid pore water saturation 
in snow with acoustic waves. The two end members correspond to no compaction 
at all and to a melting sphere packing where the grains remain in contact. 
The change of velocity for the first compressional wave was found to strongly 
depend on the compaction model, while the velocity of the second compressional 
wave mainly depends on liquid water saturation. Also the attenuation for both 
compressional waves depends mainly on the liquid water saturation and only 
little on the compaction model. Finally, a simple field experiment illustrates the 
potential use of the method to estimate the liquid water content of snow with 
acoustic sensors. 
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1. Introduction 

Knowledge of the liquid water content in snow is important to understand 
and forecast the storage and transmission rates of melt water or rain in a snow 
covered water shed m^- Also the snow albedo depends on the liquid water 
content and is important for remote sensing and energy balance modeling m- 
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This is due to the large dielectric contrast between liquid water and ice that 
increases attenuation of electromagnetic waves in the presence of liquid water. 
But the reflective behavior also depends strongly on the geometry of the pore 
space, which consists of larger grain sizes in wet or refrozen snow |39j . In fact, 
also wet snow avalanches and the corresponding stability of wet snowpacks are 
thought to depend on liquid water content, but the governing mechanisms are 
not well understood 0 EH] . 

Today, mostly electromagnetic methods are used to measure the liquid water 
content of snow in place [la m Ea ig. Surprisingly, no attempts have been 
reported to measure the liquid water content in snow with acoustic waves. 

Acoustic waves have been recognized as an important tool to characterize 
snow and acoustic wave velocities as potential index property for snow m- 
However, seemingly enigmatic observations as, for example, decreasing wave 
velocities with increasing snow density made it difficult to set acoustic observa¬ 
tions in context to other mechanical observations for snow EIIIIHI. Although 
Johnson 


and Sommerfeld m showed that Biot’s |1956a| 


1956b theory for 


wave propagation in porous materials is well suited to predict acoustic velocities 
in snow and can also explain these observations of ‘slow’ compressional waves, 
the application of acoustic waves for probing snow has remained sparse. 

A limitation for the application of Biot’s theory of wave propagation in 
porous materials is the need for a large number of properties that characterize 
the material and which are difficult to collect for snow under field conditions. 
Recently, Sidler m used empirical relationships, a priori information and ge¬ 
ometrical considerations to express snow as a Biot-type porous materials as a 
function of porosity. This model eliminates the need for a profound laboratory 
snow analysis for the application of porous wave propagation theory at the cost 
of, yet to be quantified, accuracy of the predicted results. 

Porous wave propagation theory has historically been of interest for hydro¬ 
carbon exploration and more recently also for hydrogeological applications and 
non destructive testing [24l [42l [33] . It is well known, that a mix of gas and fluid 
in the pore space significantly alters the wave propagation in a porous material 
and therefor can be used to quantify their presence [251 EH] • 

Wave phenomena tend to be sensitive to average material properties within 
a fraction of the dominant wave length. By varying the frequency of an exper¬ 
iment, acoustic waves can be used to investigate average liquid water content 
from the centimeter scale up to several meters. Therefore, the use of acoustic 
waves for the characterization of liquid water in snow has the advantage that 
the range of applicable wave lengths spans over the whole range of interest for 
spatial heterogeneity of liquid water in a snow pack. A characteristic that might 
proof to be useful to investigate the strongly heterogeneous distributions of liq¬ 
uid water commonly observed in natural snowpacks [43]. Possible applications 
may also include a combination of acoustic and electromagnetic measurements 
for higher precision without the need for calibration. 
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Table 1: Porous frame material properties for a Biot-type snow model as a function of porosity 

HQ]. 


Snow frame properties 

frame bulk modulus 

Ks (GPa) 

10 

matrix bulk modulus 

(GPa) 

30.85 

Ks{l - 4>) (7-76-0) 

shear modulus 

PS (GPa) 

3 K^{l-2v) 

2 l-\-u 

density 

p (kg/m3) 

(1-(/))• 916.7 

permeability 

K (m^) 

0 2 

'^■^(SSA)2(l-0)2 

tortuosity 

r 



2. Methods 


1956a theory of wave propaga- 


2.1. Porous material properties of snow 

For the stress strain relationships in Biot’s 
tion in porous materials four properties of a material need to be measured in 
addition to the porosity, which were named A, N, Q, and R. Biot and Willis 
[4] then introduced the poroelastic incompressibility M and the effective stress 
coefficient a to express the porous material properties in terms of the Lame 
constants A and p and the jacketed and unjacketed compressibility coefficients 
K and 5. Biot [3] then expressed the constants in terms of the bulk modulus 
of the frame material iLg, the bulk modulus of the solid matrix the shear 
modulus p, and the fluid bulk modulus Kf'm. the form which is generally used 
today. For the equations of motion the densities of the solid and fluid material 
Ps and pf, the porosity (f, and the tortuosity T of the porous material need to 
be known. The energy dissipation due to the relative motion of the fluid to the 


solid based on Darcy’s 1856 law depends on the permeability k of the material 


and the viscosity of the pore fluid. 

Therefore, seven properties define the porous frame and three define the fluid 
in the pore space. However, the individual properties are not independent of 
each other, but are based on the structure of the frame. Empirical relationships, 
a priori information, and geometrical considerations can be used to relate the 
matrix properties to each other. Snow is particularly well suited for such a 
relation because the frame matrix and pore fluids are relatively pure. Ice which 
behaves almost elastic for small strains is expected to be well suited for basic 
Biot theory where no viscoelastic effects are considered for the frame matrix 
m- For the purpose of this study we use the relations of Sidler [JD] to express 
the porous properties of snow as a function of porosity. An overview of the 
relationships that specify snow as a Biot-type porous material as a function of 
porosity are shown in Table [l] 
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Table 2: Pore fluid properties for water and air. 


Air 

density 
viscosity 
bulk modulus 

Pa (kg/m^) 
Pa (Pa s) 
ifa (Pa) 

1.29 

1.7- 10-5 
1.4- 10® 

Water 

density 
viscosity 
bulk modulus 

Pa (kg/m3) 

Pa (Pa s) 

K, (Pa) 

1000 

1.05 • 10-5 
2.25 • 10^ 


2.2. Equivalent pore fluid 

In wet snow the pore space is filled with a mix of air and water. To model 
the physics of wave propagation correctly, additional field variables would have 
to be introduced to extend Biot’s theory and account for the second phase in 
the pore space. Similar to the extension in the Double Porosity Model of Pride 
and Berryman |34L I35j . However, as a first order approximation it is common 
practice to define an effective upscaled pore fluid based on the assumption that 
one of the phases is embedded in the other phase [Mill]- As in snow there 
is only one frame matrix material and the pore space is not separated into 
two distinctly different pore regimes with individual porosity, tortuosity and 
permeability. The use of an upscaled pore fluid seems to be a valid approach 
also in this case. The effective viscosity of the pore fluid can be computed as 


?7eff = Va 



( 1 ) 


where and are the viscosity of air and water and S'w and Aa = 1 — <S'w are 
the fraction of water and air in the pore space, respectively |44] . The effective 
bulk modulus can be obtained according to the Wood equation @3 [25] 


KeS 



and the effective density is 


PeS — S'aPa + SwPw- 


( 2 ) 

( 3 ) 


The values for bulk modulus, density, and viscosity for water and air used in 
this study are given in Table 

The following steps then lead to an estimate of acoustic wave velocities and 
attenuation as a function of liquid water saturation and density of a snowpack: 

(i) Compute porosity of the snowpack from the snow density and the liquid 
pore water saturation. 
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Figure 1: The porosity in melting snow is increasing as the liquid water is accounted to 
the pore space. Porosity change as a function of pore water saturation is shown under the 
assumption of mass and volume conservation for snowpacks with different dry densities. 


(ii) Compute frame and pore fluid properties based on Table and and 
equations Q, 0, and 

(iii) Solve Biot’s equations for the obtained properties with a plane wave ap¬ 
proach (Appendix A). 


As the density of liquid water is comparable to the density of ice, the mass 
of the pore fluid can not be neglected in wet snow. The density p of the entire 
snowpack is then given by 


p= (l-(/))pi-I- </)[(!-5'w)pa + >S'wPw], (4) 

where pi, pAi and pw are the densities of ice, air and water, respectively. As 
a first order approximation, as in dry snow, the density of the air in the pore 
space (j)(l — Sv,)pA can be neglected as it is orders of magnitude smaller than 
the other densities. 

Given the density of the snow and the liquid water saturation in the pore 
space the porosity can be estimated as 

Hp,SIw) = — —(5) 

Peff Pi 

where Pes can be replaced with S„p„ as a first order approximation. The 
change of porosity with increasing water saturation for a constant density of 
the snowpack is shown in Figure 


2.3. End-member models for snow compaction 

Two end-member models for snow compaction during snow melt are pre¬ 
sented to predict how acoustic measurements would change for a fixed instal¬ 
lation in a snowpack over a period of snowmelt and to asses the need of an 


5 















independent measurement of snow density in an acoustic measurement to esti¬ 
mate the liquid water content. Existing snowpack models reproduce the com¬ 
plex interactions between meteorological conditions, geothermal heat flow and 
the corresponding response of the snow Him HI]. For this study much simpler 
models are used with the only objective to assess a realistic upper and lower 
bound for the acoustic response to an increase of liquid water content in the 
pore space. 

With the introduction of liquid water in the porous snow model shown in 
Table [T] not only the pore fluid properties change, but also the assumption that 
the density of the porous frame is also the density of the entire snowpack is 
no longer valid. Equation ([^ illustrates how this second degree of freedom is 
introduced beneath the pore water saturation. Using the end-member models 
the variation is limited to a band of plausible values. Therefore reducing the 
complexity of the analysis and leading to more comprehensive results without 
loss of generality. Moreover, if a specific property of the acoustic experiments 
shows equal variation with liquid water saturation, this is an indication that 
this property is sensitive to pore water saturation only, eliminating the need for 
an independent density measurement. 

The first end-member model assumes no snow compaction at all. It is also 
assumes that no liquid water is drained from the snowpack and therefore the 
density of the snowpack remains constant with increasing liquid water satura¬ 
tion. This end-member also corresponds to the situation were a single acoustic 
measurement on a wet snowpack is performed and evaluated against an indepen¬ 
dent density measurement. In this model the porosity increases with increasing 
water saturation 5'w following equation (§. The porosity increases because the 
part of the ice frame that has melted is accounted to the pore space. This is 
necessary as the liquid water is capable to move against the porous frame. This 
increase of porosity due to snowmelt for constant density is schematically illus¬ 
trated in Figure The striped area depicts the ice that has melted and is now 
part of the pore space. The pore space occupies now a larger fraction of the 
volume. 

There is a strong impact on wave propagation from this end-member. As the 
porosity changes, most of the properties of the porous frame change according to 
the relationships in Table In addition also the equivalent pore fluid changes 
with the increase of the liquid water saturation. The use as an end-member 
is justified by the fact, that a stronger decrease of porosity would involve an 
increase of the total volume of the snow. 

The other end-member representing a lower limit for the porosity change 
during snow melt is based on the assumption that the porosity remains constant 
due to ongoing compaction of the snowpack. This assumption can be illustrated 
with a packing of spheres where snow melt reduces the radii of the spheres and 
the snowpack compacts until the spheres are again in contact with each other. 
This end-member is illustrated in Figurej^for a cubic lattice packing of spheres. 

The constant porosity of this end-member model also leads to constant prop¬ 
erties of the porous frame. In this model, in terms of wave propagation, only 
the pore fluid properties change with increasing water saturation. While the 
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Figure 2: End-member model for constant density. Schematic illustration of the increase of 
porosity in melting snow due to accounting melted parts of the ice frame to pore space instead 
of frame. The dashed area depicts ice of the pore frame that has melted and is consequently 
accounted to the pore space due to its ability to move against the porous frame. 


porosity remains constant over the process of snow melt, the total volume of 
the snowpack reduces due to snow compaction. If the snowpack is assumed to 
compact only in vertical direction, the change of thickness of the snow cover Ah 
follows the relationship 

Ah = ho —(6) 
Po 

where po is the initial density and Ap is the change in density. The relation 
between density and liquid water saturation can be obtained using equation Q. 
For example, a dry snowpack with a density of 400 kg/m^ melting with constant 
porosity leads to a reduction of the snowpack hight of ~14 % when reaching a 
saturation of = 0.1. The use of this end-member as an upper bound of snow 
compaction is based on the consideration that a stronger increase in density 
is only possible if the geometrical constellation of the packing is changed. A 
stronger compaction is either possible if the original packing was not a dense 
packing and is changing the constellation toward a denser packing or if the 
spheres develop a dipolar distribution of radii. 

For the same amount of melted snow, the pore water saturation increases 
faster for the constant porosity end-member than for the constant density end- 
member. The percentage of the ice frame that has melted can be estimated as 


%melt 


(j)S\f/Pv/ 

(1 — 0)pi + (j)Sv/Pw 


(7) 


where the porosity ^ is a function of water saturation (j^iSw) for the end-member 
with constant density. Starting with a dry snowpack with a density of 400 kg/m^ 
{(j) — 0.56) that melts without any pore water added or removed, a pore water 
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Figure 3: End-member model for constant porosity. The illustration shows melting of a 
lattice sphere packing, a) Melting reduces the radius of the spheres, b) Snow compacts and 
ice grains stay in contact to each other. The density of the snowpack increases, but the 
porosity remains constant. 






saturation of 5'w = 0.1 corresponds to 12 % melt of the ice mass for the end- 
member with constant porosity and 16 % for the end-member with constant 
density. 

3. Results 

3.1. Change of velocity with increasing water content 

In this section phase velocities and attenuation obtained with plane wave 
solutions of Biot’s differential equations for wave propagation in porous mate¬ 
rials are shown for the presented model. The model has free degrees in pore 
water saturation, frame porosity, and frequency. These model parameters are 
inspected for an acoustic experiment to estimate liquid water content in snow 
and with consideration of the two end-member models presented above. 

A general trend of decreasing compressional wave velocity with increasing 
liquid water content in the pore space can be observed. Figure |^) shows the 
phase velocity variations at 1 kHz for the first compressional wave for snowpacks 
with different densities. Higher velocities are observed for denser snowpacks 
because of the stiffer ice frame for snow with lower porosity. Even though 
the trend is not exactly linear, the decrease of phase velocity with increasing 
water saturation is roughly 300 m/s less for an increase of 0.1 in liquid water 
saturation. This is a significant change in velocity and it should be straight 
forward to detect such velocity changes in field measurements. 

The phase velocity of the first compressional wave not only changes with 
pore water saturation but also with the density of the entire snowpack. This 
change in velocity is due to the stiffer ice frame in high density, low porosity 
snow. The increase of phase velocity is not exactly linear and varies with liquid 
water saturation as is shown in Figure]^ but is roughly an increase of 500 m/s 
for an increase of 100 kg/m^ in density. The strong gradient allows for a robust 
estimate of snow porosity in an acoustic experiment to determine liquid water 
saturation in combination with an independent density measurement. 

The attenuation for the first compressional wave is shown in Figure]^). The 
attenuation increases with increasing water content. For lighter snow the atten¬ 
uation increases stronger than for denser snow as the stiffness of the equivalent 
fluid is proportionally larger compared to the lower stiffness of the ice frame. 

Figure shows the frequency dependence of the phase velocity and atten¬ 
uation shown in Figure]^ The values for a saturation of 0.15 and a frequency 
of 1 kHz correspond to each other in these two Figures. It can be seen that 
there is only little dispersion in the phase velocities that is concentrated around 
very low frequencies. Also the attenuation is highest at low frequencies and 
exponentially decays for higher frequencies. Again, the attenuation is highest 
for lighter snow and and less pronounced for denser snow due to the stiffness 
proportions of the frame and the equivalent fluid. 

3.2. Requirement of an independent snow density measurement 

The two end-member models are used to asses the necessity of an inde¬ 
pendent density measurement to estimate liquid water content in snow with 
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Figure 4: Change of a) phase velocity of the first compressional wave at 1 kHz with increasing 
liquid water saturation in the pore space and b) the corresponding increase of attenuation. 


10 







Figure 5: Phase velocity of the first compressional wave at 1 kHz as a function of density for 
different water saturations. 


acoustic waves. An independent measurement of snow density is a limitation 
for estimating liquid water content. Automatic measurements of snow density 
are cumbersome and therefore a density measurement mostly includes manual 
work, also prohibiting self acting liquid water content measurements. In addi¬ 
tion, snowpacks can exhibit strong lateral variations in snow density. Especially 
in the presence of liquid water in the snowpack. An independent density mea¬ 
surement does not necessarily sample the same space that the remote sensing 
method does, thereby introducing an additional uncertainty. If the two end- 
member models show the same variation with liquid pore water content, such a 
measurement opens up the possibility to estimate snow density and liquid water 
content with a single acoustic experiment. 

Figure shows the decrease of compressional wave velocities as a function of 
pore water saturation for both end-member models. To compare the two end- 
members the curves are shown for a dry density of 400 kg/m^ and the density 
of the constant porosity model is adjusted with increasing liquid water content 
according to equations @ and (§. The two curves strongly spread, indicating 
that a combination with an independent density measurement is important to 
estimate the liquid water content in snow when measuring the phase velocity of 
the first compressional wave only. 

The plane wave attenuation for compressional waves of the first kind are 
shown in Figure Again the two curves for the two end-member models are 
shown to indicate minimum and maximum levels to be expected for a dry snow- 
pack density of 400 kg/m^. Especially for small amounts of liquid water in the 
pore space the two curves are close together. Indicating that the attenuation 
of the first compressional wave could be used to assess density and liquid water 
content in combination with the phase velocity. However, this is possible only if 
the geometrical effects leading to a decreasing signal amplitude are well known 
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Figure 6: Predicted a) phase velocity and b) attenuation for the first compressional wave for 
a liquid water saturation of Sw = 0.15 in the pore space. 
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Figure 7: Phase velocity of the first compressional wave for a snowpack with initial density 
of 400 kg/m® as a function of pore water saturation for the end-member model of (solid line) 
a constant porosity and (dashed line) constant density. 


and additional factors leading to attenuation as, for example, scattering of the 
signal from lateral heterogeneities can be excluded or reconstructed. 

Theory for poroelastic materials predicts a second compressional wave mode 
where the particle velocity of the pore fluid is out of phase with the particle 
motion of the solid frame. This wave mode, sometimes called the ’slow’ wave 
is difficult to observe in geological formations and is generally seen only in 
laboratory experiments [19) . Yet, the properties of snow allow, at least under 
some conditions, the propagation of this second compressional wave mode and 
the ‘slow’ wave has been measured in field conditions EH [H HO]. The second 
compressional wave is of interest because it is sensitive mostly to the properties 
of the pore space. It is therefore not surprising that in Figurej^the phase velocity 
for the second compressional wave at 1 kHz shows almost the same variation 
with increasing pore water saturation for both snow compaction end-member 
models. 

Also the attenuation of the second compressional wave that is shown in Fig¬ 
ure!^ has only little difference between the two end-member models. Interest¬ 
ingly the attenuation of the second compressional wave reduces with increasing 
liquid water content making the detection of this wave mode more prospecting. 
A potential limitation for the detection of the second compressional wave in 
wet snow is the lower porosity and accompanying higher stiffness of the frame 
that increases the attenuation of the slow wave. No observations of the second 
compressional wave are reported for densities larger than 400 kg/m^ pHl 141) . 

Due to geometrical effects and a number of mechanisms that can influence 
attenuation measurements, much more elaborate techniques are necessary to 
properly evaluate attenuation measurements compared to the measurement of 
a wave velocity. A measurement of the first and second compressional wave 
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Pore fluid water saturation (S ^) 


Figure 8: Attenuation of the first compressional wave for a snowpack with initial density of 
400 kg/m^ as a function of pore water saturation for the end-member models of (solid line) a 
constant porosity and (dashed line) constant density. 



Pore fluid water saturation {S^) 


Figure 9: Phase velocity for the second compressional wave in a snowpack with a dry density 
of 400 kg/m^ as a function of pore water saturation. The solid line corresponds to the end- 
member of constant porosity and the dashed line to the end-member of constant density. 
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Initial density n = 403.348 kg/m^ 



Figure 10: Attenuation for the second compressional wave in a snowpack with a dry density 
of 400 kg/m® as a function of pore water saturation. The solid line corresponds to the end- 
member of constant porosity and the dashed line to the end-member of constant density. 


velocity would therefor be the preferred way to estimate the density and liq¬ 
uid water content with a single acoustic experiment, supposed that the second 
compressional wave can be measured in wet snow. Numerical simulations show 
evidence that this is indeed the case under some conditions. 


4. Field example 


In this section a simple acoustic experiment in snow is used to illustrate 
the possibility to measure liquid water content in snow with acoustic waves. A 
sketch of the geometrical configuration of the experiment is shown in Figure 
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A rock was placed at the snow surface that would become an acoustic 
source later in the experiment when hit with a sledge hammer. Two snow 
pits were excavated 1 m and 3m from the location of the acoustic source and 
an accelerometer was placed 0.75 m below the snow surface in each snow pit. 
Density profiles for each snow pit were recorded by measuring the weight of a 
container with a volume of 250 cm^ that was filled by pressing it horizontally 
into the snow at 0.1 m intervals on a vertical profile. In parallel, a cylindrical 
tube with a cross-section of 30 cm^ was vertically pressed down to the bottom 
of the snow pack allowing to recover and weight a vertical snow column. The 
resulting density profiles and the average density measurements with the snow 
tube are shown in Figure 

For the experiment, the recording of the accelerometers was started with a 
sampling rate of 20 kHz and the acoustic source was excited. Figure [T^ shows 
the section of the recordings where the acoustic signal reached the accelerom¬ 
eters. The absolute zero timing was randomly chosen in front of the arrival of 
the signal at the first accelerometer. The amplitude spectra of the two recorded 
signals are shown in Figure [T^ ). The main frequency contribution of the signal 
is between 200 Hz and 800 Hz. To obtain the phase velocity and the frequency 
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Figure 11: Sketch of the acoustic experiment showing the locations of the acoustic source at 
the snow surface and the accelerometers in the two snow pits. 



Figure 12: Density profiles in the two snow pits PI and P2. The density was measured every 
10 cm and for the whole snowpack with a snow water equivalent (SWE) tube. 
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dependent attenuation the two signals were Fourier analyzed. The phase spec¬ 
trum was used to compute the phase velocity and the amplitude spectrum was 
used to compute the attenuation [37l [301E21ES] . The corresponding results are 
shown in Figure [l4^) and 14 d), respectively. 

As there is little to no dispersion expected in this experiment (see Figure 
the phase velocity is best evaluated at the central frequency of the propa¬ 
gating pulse. Which is around 500 Hz for this experiment. The phase velocity at 
other frequencies are used to obtain an estimate on the precision of the velocity 
measurement. Based on Figure [l^ ) the maximum and minimum velocity for 
wave propagation between the two sensors was chosen to be 1460 m/s and 1570 
m/s, respectively. 

To obtain an estimate for the volume fraction of liquid water content the 
phase velocity and the independent density measurement were inverted for 
porosity </> and liquid water saturation of the pore space S^i as follows: 


1. Use a least square cost function equally weighting errors in velocity and 
density and chose random initial values for the unknown porosity and 
liquid water saturation. 


2. Calculate the properties of a Biot-type porous material and an equivalent 
pore fluid for the chosen initial values using Table and and equations 
0, 0, and 0. 


3. Compute the phase velocity for the resulting material for the plane wave 


solution of Biot’s equation for wave propagation in porous materials (Ap¬ 


pendix A). 


4. Evaluate the cost function and optimize the initial values by repeating 
steps 2 to 4 for adjacent values of the initial guess according to the Nelder- 
Mead scheme [5S]. 


The resulting porosity and liquid water saturation as well as the correspond¬ 
ing volume fraction of liquid water content are shown in Table Depending 
on the velocity and density measurements the liquid water content for the ex¬ 
amined snowpack is estimated to be between 4 % and 8 % and can be classified 
as wet snow (3) on a scale from (1) which corresponds to 0 % Vol to (5) which 
corresponds to more than 15 % Vol [13]. This corresponds with the observa¬ 
tions that liquid water was discharging from the snowpack, but water could 
not be pressed out of the snow by squeezing the snow in by hands. No alter¬ 
native quantitative measurements of liquid water content have been made and 
the rule of thumb observations are an admittedly weak source for comparison. 
However, it is interesting to see that the estimate for the porosity is relatively 
constant in the range of the estimated precision of the measurements and that 
the liquid water content depends more strongly on the precision of the density 
measurement than on the phase velocity estimate. 

In future research the liquid water estimates with acoustic methods will be 
compared to quantitative measurements with alternative methods. Also the spa- 
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Figure 13: Normalized a) acceleration and b) corresponding amplitude spectra for the mea¬ 
surements in the two snow pits. The acceleration of the second receiver is corrected for 
spherical spreading and multiplied by 1.7 to have the same maximal amplitude as the first 
receiver. 
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Figure 14: Estimated a) phase velocity and b) attenuation from the measurements shown 
in Figure The amplitudes were corrected for spherical geometrical spreading to compute 
the attenuation. The minimum and maximum values evaluated in Tableware indicated with 
horizontal lines. 
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Table 3: Estimated snow porosity and liquid water content for minimum and maximum 
estimated phase velocities and densities. 



Vest = 1460 


Vest = 1570 


p = 408 kg/m^ 

^ = 0.50 

5 % Vol 

= 0.10 

= 0.48 

4 % Vol 

= 0.07 

p = 552 kg/m^ 

(p = 0.50 Sw 
9 % Vol 

= 0.18 

(j) = 0.48 

8 % Vol 

= 0.15 


tial sensitivity of the measurement and the strongly heterogeneous distribution 
of the liquid water in the pore space will be investigated. 

5. Conclusions 

An equivalent pore fluid in combination with empirical relationships for a 
porous snow frame model are used to asses the plane wave solutions of Biot’s 
theory for liquid water in snow. The results indicate that acoustic waves are 
sensitive to liquid water content in snow and the phase velocity of the first 
compressional wave is expected to be lowered by roughly 300 m/s for every 
increase of 0.1 in water saturation in the pore space. In combination with an 
independent density measurement acoustic waves should be capable to estimate 
the liquid water content in snow. The use of two end-member models for snow 
compaction reveal that additional measurements of attenuation of the first com¬ 
pressional wave or the phase velocity of the second compressional wave, should 
this be possible, would allow to estimate also the density of the snowpack from 
the acoustic measurements. A field measurement illustrates the simplicity and 
robustness of the method and highlights the importance of the precision of the 
density measurement on the estimate of liquid pore water content. 
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Appendix A. Compressional plane wave velocities in a Biot-type 
porous material 


A convenient method to solve Biot’s 1956b differential equations for wave 
propagation in porous materials in closed form is to assume plane waves as a 
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(A.l) 

(A.2) 


solution nig. Plane waves of the form 

^s = ^so 

^f='dfo 

where i?s and are the velocities of the waves propagating in the solid and fluid, 
to is the angular frequency, k is the wavenumber and i = V—1, are substituted 
into Biot’s differential equations 

= PiidtA^^ + Pi2dttuA + KuA - (A-3) 

-<t>djPf = Pi2dttuA + P22dluA - b{uA - uA), (A.4) 

where a is the stress tensor, pij are mass coefficients which take into account that 
the relative fluid flow through the pores is not uniform, b = (j)^r]/K is the friction 
coefficient with rj being the viscosity, (/> the porosity, and n the permeability, u is 
the displacement vector with the superscripts (S') and (/) referring to the solid 
matrix and the fluid, respectively. Equations |A.3| and |A.4| have to be expressed 
in the velocity stress formulation to obtain the quadratic dispersion relation 


-[Pj + -Yp 


-Y Kg + -ps + M{2apf - p) 


M + -psj = 0, (A.5) 

that can be solved for the two complex plain wave moduli Vc = uj/ k of the first 
and second compressional wave. Symbols in equation A.5 are the fluid density 
Pf, the viscodynamic operator Y 


Y{t) = mdtS{t) + —6{t) 


with S{t) being the Dirac function and 


m = 


PfT 


where T is the tortuosity. Gassman’s 1951 modulus Kq is 

Kg = K^+ a^M, 
with the effective stress coefficient 

Km 


a = 1 — 


Ks 


the poroelastic incompressibility 

M = 


(« - 4>) 

Ks Kf 


-[ -1 


(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.IO) 
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and Ks, Km and Kf being the frame, matrix and fluid bulk moduli. The phase 
velocity and attenuation can be obtained from the complex plain wave modulus 
using 


= [Re(KM)-i] 


(A.ll) 


and 



Qp{u}) ^ 


g lm(V'cM) 

Re(K(a;))’ 


(A.12) 
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